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Prospects for indirect detection of sneutrino dark matter with IceCube 
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We investigate the prospects for indirect detection of right-handed sneutrino dark matter at the 
IceCube neutrino telescope in a U{\)b-l extension of the MSSM. The capture and annihilation of 
sneutrinos inside the Sun reach equilibrium, and the flux of produced neutrinos is governed by the 
sneutrino-proton elastic scattering cross section, which has an upper bound of 8 x 10 -9 pb from 
the Z' mass limits in the B — L model. Despite the absence of any spin-dependent contribution, 
the muon event rates predicted by this model can be detected at IceCube since sneutrinos mainly 
annihilate into leptonic final states by virtue of the fermion B — L charges. These subsequently 
decay to neutrinos with 100% efficiency. The Earth muon event rates are too small to be detected 
for the standard halo model irrespective of an enhanced sneutrino annihilation cross section that 
can explain the recent PAMELA data. For modified velocity distributions, the Earth muon events 
increase substantially and can be greater than the IceCube detection threshold of 12 events km -2 
yr _1 . However, this only leads to a mild increase of about 30% for the Sun muon events. The 
number of muon events from the Sun can be as large as roughly 100 events km -2 yr -1 for this 
model. 



PACS numbers: 12.60.Jv, 95.35.+d, 14.60.Lm 



I. INTRODUCTION 



There are various lines of evidence supporting the ex- 
istence of dark matter in the universe, but its identity 
remains a major problem the solution to which likely 
rests at the interface of particle physics and cosmology. 
It is well established that particle physics can explain 
dark matter in the form of weakly interacting massive 
particles (WIMPs) [1[ . In the standard scenario the dark 
matter relic abundance, as precisely measured by cosmic 
microwave background (CMB) experiments Q is deter- 
mined from the thermal freeze out of dark matter anni- 
hilation in the early universe. There are currently major 
experimental efforts for direct and indirect detection of 
dark matter particles. Indirect detection investigates an- 
nihilation of dark matter to various final states (photons, 
anti-particles, neutrinos) through astrophysical observa- 
tions, while direct detection probes the scattering of the 
dark matter particle off nuclei inside dark matter detec- 
tors. 

Supersymmetry is a front-runner candidate to address 
the hierarchy problem of the standard model (SM). The 
minimal supersymmetric standard model (MSSM) has 
become the focus of major theoretical and experimental 
activities for the past two decades. It has a natural dark 
matter candidate, namely the lightest supersymmetric 
particle (LSP), which can have the correct thermal relic 
abundance ||. It is also believed that there are gauge 
symmetries beyond those of the SM. A minimal exten- 
sion of the SM gauge group, motivated by the nonzero 
neutrino masses, includes a gauged U(\)b-l gauge sym- 
metry [i[ (B and L are baryon and lepton number respec- 
tively). Anomaly cancellation then implies the existence 
of three right-handed (RH) neutrinos and allows us to 
write the Dirac and Majorana mass terms for the neutri- 



nos to explain the light neutrino masses and mixings. 

The B — L extended MSSM also provides new dark 
matter candidates: the lightest neutralino in the B — L 
sector Q and the lightest RH sneutrino In this 
work we will focus on the sneutrino as the dark matter 
candidate 1 . The candidate is made stable by invoking a 
discrete i?-parity, but in the context of a B — L symmetry, 
a discrete matter parity can arise once the U{\)b-l is 
spontaneously broken [9J] . The B — L gauge interactions 
can yield the correct relic abundance of sneutrinos if the 
U(1)b-l is broken around the TeV scale. 

Recently, it has been shown that it is possible to 
explain the positron excess observed in the PAMELA 
data [lCj in the context of a low scale B — L extension 
of the MSSM @, [n, El- Due to a factor of 3 difference 
between the B — L charges of the quarks and leptons, the 
anti-proton flux is naturally suppressed in this model in 
agreement with the PAMELA anti-proton data. Further- 
more, the U(1)b-l gauge coupling unifies with those of 
the SM symmetries, and the B — L symmetry can be bro- 
ken radiatively. The B — L breaking around a TeV results 
in a Z' gauge boson with around a TeV mass that can 
be probed at the LHC along with the other new states of 
this model. 

The RH sneutrino of this B — L extended model can 
be detected when it elastically scatters off a nucleus. The 
sneutrino-proton scattering cross section is large enough 
to be probed in the ongoing and upcoming dark matter 



1 It is also possible to have successful inflation in the context of 
the U(1)b-l model In this case the dark matter candidate 
(the RH sneutrino) can become a part of the inflaton field and 
thereby gives rise to a unified picture of dark matter, inflation 
and the origin of neutrino masses 0. 
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direct detection experiments Q. In addition, annihila- 
tion of sneutrinos at the present time produces LH neu- 
trinos. It is interesting to investigate the possibility of 
indirect detection of sneutrino dark matter by using final 
state neutrinos in the IceCube neutrino telescope. This 
ongoing experiment plans to probe the neutrino flux aris- 
ing from the annihilation of gravitationally trapped dark 
matter particles in the Sun and the Earth. We will ex- 
amine the status of the U(1)b-l model in two cases. In 
case 1, the sneutrinos annihilate mostly into RH neutri- 
nos that subsequently decay into LH neutrinos and the 
MSSM Higgs. In case 2, the sneutrinos annihilate mostly 
into the lightest Higgs boson in the B — L sector that 
decays into t + t~ and b + b~ quarks, which subsequently 
produce LH neutrinos via three-body decays. The recent 
PAMELA data [l(| can be explained in case 2, where 
the final state taus give rise to the positron excess in 
the cosmic ray flux without producing a significant num- 
ber of antiprotons @, [ill, E3]- The large cross section 
required to explain the data arises from Sommerfeld en- 
hancement [Lf or from the non-thermal production of 
dark matter [12|. 

Since the source of neutrinos are different in the two 
cases, two-body versus three-body decay, the energy 
spectrum of the neutrinos can be used to distinguish the 
cases. We will estimate the muon neutrino flux as well as 
the muon flux in both scenarios as a function of sneutrino 
mass. Since the Large Hadron Collider (LHC) is on the 
verge of producing physics results, it will enable us to 
measure the mass of the dark matter candidate. There- 
fore, using the LHC measurements and the IceCube re- 
sults in tandem, we hope to discern the B~L model. We 
will present predictions of this model using the standard 
dark matter halo model as well as the modified velocity 
distributions obtained in recent galaxy simulations. 

This paper is organized as follows. In section II, we 
discuss the low scale U(1)b-l model. In section III, we 
give a general discussion of the indirect detection of sneu- 
trino dark matter via neutrino final states. In section IV, 
we present our results and discuss the prospect of detec- 
tion of sneutrino dark matter at IceCube in case 1 and 
case 2. In section V, we show the results obtained for 
the modified velocity distributions. In section VI, we 
compare predictions for the sneutrino dark matter in the 
U(1)b-l model with those for the neutralino dark mat- 
ter in the minimal supergravity model. Finally, we close 
by concluding in section VII. 



II. THE U{1) B -l MODEL 

Since this B — L is a local gauge symmetry, we have 
a new gauge boson Z' (and its supersymmetric partner) . 
In the minimal model, we also have two new Higgs fields 
H[ and H 2 (that are SM singlets) and their supersym- 
metric partners. The vacuum expectation values (VEVs) 
of these Higgs fields break the B — L symmetry. We 
can write the superpotential of the model as follows (the 



boldface characters denote superfields) 

W = Wmssm + Wb-l + 2/dNH u L , (1) 

where H u and L are the superfields containing the Higgs 
field that gives mass to up-type quarks and the LH lep- 
tons respectively. For simplicity, we have omitted the 
family indices. The Wb-l term contains H^, H' 2 and 



N [llj. Its detailed form depends on the charge assign- 
ments of the new Higgs fields. 

The U{\) B -l is broken by the VEV of H[ and H 2 , 
which we denote by v[ and v' 2 respectively. This results in 
a mass mz' for the Z' gauge boson. We have three phys- 
ical Higgs fields <j), $ (scalars) and A (a pseudo scalar). 
The masses of the Higgs fields follow mi < cos 2 (2/3')m|/ 
(where tan/3' = (H 2 ) / (H[)) and to$, m.4 ~ mz>. 

Various B — L charge assignments are allowed by 
anomaly cancelation. We choose the charge assignment 
shown in Table 1 . In this case H 2 couples to the RH neu- 
trinos and gives rise to a Majorana mass upon sponta- 
neous breakdown of the U (1)b-l- Choosing these Majo- 
rana masses in the 100 GeV — 1 TeV range, we have three 
(dominantly RH) heavy neutrinos and three (dominantly 
LH) light neutrinos. The masses of the light neutrinos 
are obtained via the see-saw mechanism. 



Fields 


Q 


Q c 


L 


L c 


H[ 




Qb-l 


1/6 


-1/6 


-1/2 


1/2 


1 


-i 



TABLE I: The B — L charges of the fields for the minimal 
model. Here Q and L represent quarks and leptons respec- 
tively, while H[ and are the two new Higgs fields. The 
MSSM Higgs fields have zero B — L charges. 

A natural dark matter candidate in this model is the 
lightest sneutrino N . We note that it has fewer gauge 
interactions than other supersymmetric particles, and its 
mass receives the smallest contribution from the gaugino 
loops. Based on the dominant channel for sneutrino an- 
nihilation we therefore consider the following two cases: 

• Case 1: A generic case where a solution to 
the positron excess observed by PAMELA is not 
sought. In this case the dominant annihilation 
channels are the S'-wave processes NN — > NN 
and N*N* — > N*N* via f-channel exchange of 
Z'. There are also NN* -► NN*, ff annihila- 
tion modes via s-channel exchange of a Z' or B — L 
Higgs fields, but these are P-wave suppressed and 
can be completely neglected (particularly at the 
present time). In this case the annihilation cross- 
section has the nominal value ~3x 10 -26 cm 3 /sec 
(dictated by thermal freeze out) at all times. The 
RH neutrinos produced from dark matter annihila- 
tion quickly decay to LH neutrinos and the MSSM 
Higgs. 

• Case 2: In this case the PAMELA puzzle is ad- 
dressed via Sommerfeld enhancement of sneutrino 
annihilation at the present time [11| . In this part of 
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the model parameter space the lightest B — L Higgs 
<f> is much lighter than the Z' . The dominant an- 
nihilation channel is N*N — > cjxfi via the s-channel 
exchange of the <fi or <E>, the t or u-channel exchange 
of a N, and the contact term \N\ 2 tfi 2 . The inter- 
actions for these processes arise from the £>-term 
part of the potential, and their strength is propor- 
tional to my. There are other S- wave processes 
with Higgs final states N*N -> 0$, <j)A, A4, 
but they are kinematically suppressed and/or for- 
bidden. The annihilation modes NN —* NN and 
N*N* — > N*N* are also subdominant in this case. 
As in the previous case, annihilations to // final 
states are P-wave suppressed and hence totally neg- 
ligible. The cross section for annihilation to the 
4>4> final state at the present time is required to be 
3 x 10 -23 cm 3 /sec in order to explain the PAMELA 
data. Sufficient Sommerfeld enhancement is ob- 
tained as a result of the attractive force between 
sneutrinos due to the 4> exchange provided that 
the mass of (p is small (< 20 GeV) 2 . The tj> subse- 
quently decays into fermion-antifermion pairs very 
quickly via a one-loop diagram, and it mostly pro- 
duces t + t~ final states by virtue of the fermion 
B — L charges fill ]. 

The sneutrino-proton scattering cross section for this 
model can be in the 10~ n — 10~ 8 pb range for a reason- 
able choice of parameters that satisfy the relic density 
constraint, cf. 0, [HI- This opens up the prospect for 
direct detection with the help of the next generation of 
experiments 14| . The current upper bound for the spin- 
independent cross section is 4.6 x 10 ~ 8 — 2 x 10 ~ 7 pb for 
a dark matter mass of 60 — 1200 GeV; this is just above 
the highest possible values for our model 3 . 



boson 4 . We assume for most of this paper that the to- 
tal LH neutrino flux branches into every neutrino flavor 
equally (see subsection IIV Al for a discussion) . Assum- 
ing that the mass difference between the RH sneutrinos 
and RH neutrinos is small 5 , the RH neutrinos are pro- 
duced non-relativistically, and hence each LH neutrino 
and Higgs receives an energy equal to half of the sneu- 
trino mass. 

In case 2, RH neutrinos constitute about 10% of the 
annihilation final states. Two of the lightest B — L Hig- 
gses 4> compose the remaining 90% of the branching frac- 
tion. This branching fraction is necessary to provide a 
high enough leptonic particle rate to fit the PAMELA 
data. As mentioned in the previous section we need 
< 20 GeV. For 4 GeV < m < 20 GeV, the final 
states are mostly taus (74%) and b quarks (16%), where 
the dominance of tau final states is a result of the fermion 
B — L charges. The LH neutrinos in this case arise from 
the three-body decay of taus and bottom quarks. For 
77i0 < 4 GeV, we would have mostly muons and charm 
quarks. 

Both the case 1 and case 2 scenarios of our model 
display a crucial signature difference when compared to 
the standard neutralino LSP in the MSSM. The energy 
distribution of the produced LH neutrinos from the RH 
neutrino decay is a delta function occurring at half of 
the sneutrino mass. Other annihilation channels in this 
model, as well as those available in the MSSM, produce 
additional neutrino signal via three-body decays such as 
t~ — > e~v T v e . This difference opens up a significant pos- 
sibility to differentiate between the B — L model and the 
MSSM with the help of the differential energy spectrum 
of the detector event rates. This is discussed further in 
section IIVI 



B. Neutrino Flux 



III. PROSPECTS FOR INDIRECT DETECTION 
AT ICECUBE 

A. The Neutrino Signal 

The B — L model also shows great promise for indi- 
rect detection, and we focus in particular on the potential 
neutrino signal at the IceCube experiment. In case 1, the 
sneutrinos annihilate to produce RH neutrinos that sub- 
sequently decay into a LH neutrino and a neutral Higgs 



Sneutrino annihilation in the Sun and the Earth pro- 
duces an expected neutrino flux through IceCube. This 
flux is modeled by calculating the number of gravitation- 
ally captured sneutrinos and then considering the prop- 
agation and detection of the produced neutrinos. The 
number of captured dark matter particles as a function 
of time is governed by a differential equation the solution 
to which is 



N(t) = \j — tanh VCAt 



(2) 



2 It is possible to invoke a non-thermal scenario where the sneutri- 
nos are created from the decay of heavy moduli or gravitinos |12H . 
In this case we do not need Sommerfeld enhancement to satisfy 
the PAMELA data, and the annihilation cross section will be 
large, 3 X 10 — 23 cm 3 /sec, at all times. 

3 Since B — L symmetry is vectorial, the spin-dependent cross sec- 
tion is zero in this model. 



4 RH neutrino decay to a charged lepton and a charged Higgs is 
typically forbidden. 

J This is the case when the soft supersymmetry breaking mass 
of the sneutrino is similar to or smaller than supersymmetry 
conserving Majorana mass of the (s)neutrino. A rather small 
soft mass term is motivated if the B — L symmetry is to break 
radiatively and is needed to keep the lightest B — L Higgs <j> light 
as in case 2 [Til . 
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where C is the total capture rate and depends on both 
the total scattering cross sections off nucleons and A is 
related to the annihilation cross section; see Ref. [17| for 
details. The total rate of annihilation is given by 



C 

Ta = 77 tanh 



t 



1 eq 



(3) 



The number of captured sneutrinos will saturate as long 
as the length of time for t he p rocess has exceeded the 
equilibration time, r eq = (V CA)~ l . 

In equilibrated systems, the rate of annihilation is en- 
tirely dominated by the capture rate C, Ta ~ C/2. We 
can explain equilibration in the B — L model by consid- 
ering some example cross sections. Since the age of the 
solar system is 4.5 Gyr, for a 1 TeV sneutrino with an 
annihilation cross section of 3 x 10 -23 cm 3 /sec (3 x 10~ 26 
cm 3 /s), a spin- independent cross section ctsi of at least 
10~ n pb (10~ 8 pb) is needed to reach equilibration in 
the Sun. This assumes no spin-dependence as the B — L 
model has none. The scattering cross section needed to 
achieve equilibration in the Earth is already excluded by 
direct detection bounds. 

Alternatively we can fix the scattering cross section 
and place a limit on the annihilation cross section. In 
the B — L model, the cross section for sneutrino-proton 
elastic scattering follows 



a S i oc 



( 9b-lQl 



(4) 



where gs-L and Ql are the U(1)b-l gauge coupling 
and B — L charge of leptons, respectively, and m p is the 
proton mass. The limits on the Z' mass from LEP and 
Tevatron are given by (TBI . HE] , 



mz> 



9b-lQl 



> 6 TeV 



(5) 



This results in an upper limit on crgi of 8 x 10~ 9 pb. 
Assuming this bound is realized, an annihilation cross 
section > 4 x 10~ 26 cm 3 /s (1 x 10~ 18 cm 3 /s) needs to be 
achieved to reach equilibrium in the Sun (Earth). Note 
that we can always choose the B — L gauge coupling and 
scale B — L charges in accordance with anomaly cancela- 
tion such that asi is saturated while obtaining the correct 
relic density for sneutrino dark matter. This is possible 
since a different combination of gs-L and Qb-l appears 
in the relic density calculation . This is in contrast to the 
MSSM case where the SM gauge couplings and charges 
are fixed. 

Since equilibrium is easily achieved in the Sun, the 
neutrino signal will depend solely on C, or equivalcntly 
ctsi, so the increased annihilation rate in case 2 of our 
model confers no advantage compared to typical MSSM 
cases for annihilation in the Sun. On the other hand, 
choosing reasonable values for either of the relevant cross 
sections demonstrates that equilibrium is nearly impos- 
sible to reach for the Earth without significant deviation 



from the assumptions made in [17]. Consequently, the 
neutrino signal from the Earth will depend on both C 
and A. Therefore one expects a much larger signal for 
case 2 as compared to either case 1 or the neutralino 
dark matter models [H| . 

The annihilation of sneutrinos in the Sun and Earth 
yields neutrinos that can be detected by the IceCube 
experiment. IceCube can distinguish between neutrino 
signals from the Earth and Sun with the help of an angle 
cut. This cut restricts the detection to an angle range 
of 90° < 9 < 113° in the case of the Sun, where 9 is 
the Earth zenith angle. One has to measure below the 
horizon to be able to distinguish the background of at- 
mospheric neutrinos from the signal, and the Sun cannot 
be more than 23.5° below the horizon at the South Pole 
[27l |28| . In the case of a search for a potential Earth sig- 
nal one looks at a zenith angle of about 180°, i.e., directly 
to the core of the Earth [271 ] . 

Muon neutrinos create muons via charged current in- 
teractions in the detector. The qualitative behavior of 
the muon flux depends on the corresponding neutrino 
muon flux, and the differential neutrino spectrum is given 
by 



dN v 

dE v AttD 



r>2 n 



N dE,, 



(6) 



see for example Ref. [2l|. AppendixlA~lcontains a detailed 
discussion about the mass dependence of this equation. 

The IceCube detector records the Cerenkov light from 
relativistic charged particles in its volume. Cosmic ray 
showers create a muon background signal that can be 
controlled by selecting for upward-going and contained 
muon events. The atmospheric neutrino background is 
well understood and may be subtracted away from the 
signal. 

In addition to the muon flux through IceCube, elec- 
tromagnetic and hadronic cascades inside the detector 
might also allow sneutrino dark matter detection. Elec- 
tromagnetic cascades occur via charged current inter- 
actions. By depositing some of the incoming neutrino 
energy in taus and electrons, Bremsstrahlung radiation 
produces a localized cascade of energy that the digital 
optical modules of IceCube can record. In the results of 
Appendix [B] we have ignored any contribution from the 
charged current electromagnetic cascades of the muons, 
since their contribution has already been considered in 
the form of Cerenkov radiation from the muon tracks. 
Hadronic cascades occur for both neutral current and 
charged current interactions. As the neutrino scatters 
off of a nucleus in the detector, the nucleus breaks up 
and produces products such as pions that in turn decay 
into detectable photons. Note that for neutral current 
interactions the energy of the outgoing neutrino is lost 
and is not recorded in any cascade. The energy from 
localized electromagnetic and hadronic cascades is much 
harder to reconstruct compared to muon tracks but still 
might produce an interesting signal in the detector, see 
Appendix [B] 
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IV. MODEL RESULTS 

The annihilation of sneutrinos in the Sun and Earth 
results in a flux of particle events through the IceCube 
detector that are calculated using DarkSUSY, which uses 
results from WimpSim [U [13] • The calculations account 
for neutrinos produced via decays, as well as neutrino os- 
cillation, loss via charged current interactions and scat- 
tering via neutral current interactions. DarkSUSY de- 
fault parameters are used, which include a Gaussian dark 
matter velocity distribution and an NFW halo profile. 
Realistic Sun and Earth density profiles are integrated 
over numerically according to [231 ] . For both case 1 and 
case 2, the maximum spin-independent cross section al- 
lowed by the Z' limits is used. Similarly, the annihilation 
cross section is fixed at 3 x 10~ 26 cm 3 /s (3 x 10~ 23 cm 3 /s) 
for case 1 (case 2). Finally, the results presented in the 
subsections below use the convention of a detector energy 
threshold of 1 GeV. IceCube effective areas have not been 
calculated for our model, but we anticipate that they 
would be slightly larger than those used for the MSSM 
scenarios since we have a slightly harder spectrum. This 
is especially true in case 1. 



A. Sensitivity to Neutrino Flavor 

For the results that follow we have considered equal 
branching to the three flavors of LH neutrinos, but in 
principle this need not be the case. The exact flavor 
composition of LH neutrinos produced from sneutrino 
annihilation in the Sun depends on the detailed structure 
of Majorana and Dirac couplings in the neutrino sector. 
In Fig. [TJ the resulting muon neutrino flux for a 100% 
branching ratio to a single flavor is compared to equal 
flavor ratios in both case 1 and case 2 (upper and lower 
panels respectively). 

It is seen from the upper panel that in case 1 for 
sneutrino masses below 300 GeV (LH neutrino energy 
below 150 GeV) flavor composition of produced neutri- 
nos does not matter since oscillations are very efficient at 
low energies and easily mix the neutrino flavors. There- 
fore 100% v e , is^, or v T each leads to the same signal at 
the detector. However at high energies oscillation length 
Lose °c E„/Am 2 elongates, and oscillations become less 
efficient. Here Am 2 is the difference between (mass) 2 
of neutrino mass eigenstates. This effect is most impor- 
tant for v e 's since they oscillate to u^'s via the small mass 
splitting responsible for solar neutrino oscillations Am 2 ol . 
This is why the flux at the detector falls quickly for 
100% v e branching ratio at high energies. The effect is 
less pronounced for 100% and v T branching ratios be- 
cause the relevant mass splitting is the one responsible for 
atmospheric neutrino oscillations Am 2 tm , which is much 
larger. However, it is seen that the flux for 100% 
Vp branching ratio is less than that for 100% v T branch- 
ing ratio at high energies. This is because of charged 
current interactions inside the Sun whose cross section 
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(b) Sun Muon Neutrino Flux, Case 2 

FIG. 1: Total muon neutrino rates received at the Earth for 
the U(1)b-l model as a function of the sneutrino mass in 
the case of sneutrino dark matter capture and annihilation in 
the Sun. The results are for one year of detection with Ice- 
Cube. The B — L model is robust to changes in the neutrino 
branching ratios. 100% branching to i/ e , and v T is shown 
in orange (bottom line in case 1), green (second line from 
the bottom in case 1) and blue (top line in case 1) respec- 
tively. Results of equal branching to neutrino flavors are in 
red (second line from the top in case 1). 



is proportional to the neutrino energy. These interac- 
tions convert muon neutrinos to muons that are quickly 
stopped in the Sun due to electromagnetic interactions 
that result in attenuation of the neutrino flux. Charged 
current interactions also convert tau neutrinos to taus. 
However, due to their much shorter lifetime, they decay 
back to v T before any significant energy loss. Neverthe- 
less, for sneutrino masses up to 1.5 TeV, the result for 
equal branching ratios to three flavors is within a factor 
of a few compared with the 100% branching ratio to a 
single neutrino flavor. Moreover, for a typical model, it 
is unlikely that sneutrino annihilation produces only one 
flavor of RH neutrinos. Therefore equal branching to the 
three flavors is a good approximation in case 1. 

In case 2, the lower panel 6 , there is virtually no differ- 



6 The effect of the 1 GeV conventional energy threshold in the 
spectrum can be seen at low masses as more of the neutrino signal 
is lost under the threshold; this causes the maximum event rate 
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ence between various flavor compositions. This is because 
sneutrino annihilation mainly produces taus in this case 
(the branching ratio for production of RH neutrinos is 
only 10%). Hence equal branching to the three flavors is 
a nearly perfect approximation in this case. 

We conclude that our results do not depend critically 
on the choice of neutrino flavor branching ratios in cither 
case. 



B. Contributions to Muon Flux 

It is worth emphasizing that case 1 and case 2 
yield different neutrino signals. In case 1, LH neutri- 
nos are produced from two-body decay of (almost non- 
relativistic) RH neutrinos. This produces a delta func- 
tion in the energy of the LH neutrinos at one-half the 
mass of the sneutrino dark matter 7 . On the other hand, 
in case 2, the sneutrino dominantly annihilates to <jxf> 
final states, and each <jj decays to a fermion-antifermion 
pair via a one-loop diagram. The partial decay rate of <p 
is proportional to the squares of the mass of the resulting 
fermion and the fourth power of its B — L charge 0, [ll[ . 
As a result, the largest contribution to the annihilation 
is from taus (ss 74%) and bottom quarks (« 16%), where 
the quark signal is suppressed due to the B — L charge. 
Both of these final states produce neutrinos via three- 
body decay that results in a spread in energy signal. 



1.5x10 



Fig. 2(a) shows the muon neutrino flux energy spec- 



trum through a kilometer squared of IceCube in one year 
for a 300 GeV sneutrino for case 1. The delta function 
at half the mass of the sneutrino can be seen clearly. A 
small portion of muon neutrinos from this initial anni- 
hilation state are scattered via neutral current interac- 
tions inside the Sun to lower energies. This produce s the 
slight bump in the spectrum at low energies. Fig. |2(b)| 
plots the resulting muon flux from the charged current 
interactions inside the IceCube detector. As expected 
for a monochromatic incident neutrino, the spectrum of 
muons has a linear dependence on energy. 

For case 2, the delta function from the neutrino chan- 
nel at the detector is subdominant to the other annihi- 
lation channels, see Fig. 3(a) First, consider that the 
sneutrino annihilation mainly produces taus and bot- 
tom quarks that subsequently produce LH neutrinos via 
three-body decays. Second, due to the larger sneutrino 
mass of 1 TeV (in order to explain the PAMELA data), 
the LH neutrinos produced from two-body decays have 
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FIG. 2: In the upper (lower) panel, muon neutrino (muon) 
flux through IceCube from annihilation of 300 GeV sneutrinos 
in the Sun for case 1. 



a higher energy than in the 300 GeV case. Therefore 
they lose energy via neutral current interactions and get 
absorbed via charged current interactions inside the Sun 
more efficiently. As a result of both of these facts, there 
are more neutrinos with low energies at the detector from 
each channel in this case than in case 1. This also is re- 
flected in the spectrum of muon flux, shown in Fig. 3(b)| 
which does not show a linear dependence on energy due 
to the presence of three-body decays. This is in contrast 
to Fig. 



2(b) 



The muon event signal from annihilation in the Earth 
for case 1 and case 2 is too small to detect since the dark 
matter population has not reached equilibrium; there- 
fore, the production of neutrinos depends on both the 
scattering cross section and the annihilation cross sec- 
tion, which is small in this scenario. 



to move to the right from the edge of the graph. This effect is not 
evident in case 1 since the majority of the neutrino flux arrives 
at higher energies and is unaffected by the small threshold. 
7 There is one additional potential source for neutrinos: the Higgs 
produced from the decay of the RH neutrinos can itself decay 
to a bb pair. We checked that this contribution gives only a 
few percent change in the signal. We therefore neglect it in our 
numerical calculation for the sake of simplicity. 



C. Mass Dependence of Muon Flux 

Fig. 2] shows our results for the total muon rate inte- 
grated over energy as a function of the sneutrino mass 
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FIG. 3: The same as Fig. [2] but with a 1 TeV sneutrino in 
case 2. Individual annihilation channels are shown: neutrino 
(red, dotted), tau (green, dashed), bottom quark (purple, dot- 
dashed) and all channels (blue, solid). 

rrijy for annihilation in the Sun 8 . The figure shows 
both the case 1 and case 2 rates in events km~ 2 yr _1 . 
The plots have two characteristics: an increase at lower 
masses culminating in a peak followed by a general de- 
crease in event rates at higher masses. 

The decrease of the event rates for higher is re- 
flective of the decrease of the neutrino flux due to the 
kinematic suppression of sneutrino capture (the factor 
scales approximately like l/m^ for large masses 9 ). The 
linear increase at low is explained by the linear de- 
pendence of the cross section for charged current inter- 
actions on the energy of neutrinos at the detector (which 
is proportional to the sneutrino mass). The case 1 sig- 
nal is larger than the case 2 signal for lower values of 
sneutrino mass. LH neutrinos are produced in two-body 
decays in case 1 versus three-body decays in case 2, and 
hence have a higher energy. As a result, the cross section 
for conversion of neutrinos to muons at the detector is 



8 The apparent discrete nature of these plots occurs because only 
a few values of sneutrino mass are recorded in the WimpSim ta- 
bles used by DarkSUSY; the program interpolates between these 
points. The effect is numerical and not physical. 

9 See Appendix lAl for a more detailed definition of "large". 
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FIG. 4: Total muon rates detected at the Earth from 
annihilation of sneutrino dark matter in the Sun as a 
function of the sneutrino mass. The results are for one 
year of detection with IceCube. Case 1 (case 2) is the 
highest (lowest) peaked line. The dotted line denotes the 
mass range where one cannot explain the PAMELA data 
using case 2 anymore. 



larger in case 1. However, for large sneutrino masses 
case 1 has a smaller signal than case 2. The produced 
LH neutrinos, 100% of case 1 products, get absorbed 
via charged current interactions or lose energy via neu- 
tral current interactions inside the Sun more efficiently 
because of their larger energy, thus a smaller number of 
neutrinos arrive at the detector. 

Refs. [lH, display sensitivity plots for the detec- 
tion of a muon signal in the case of standard neutralino 
dark matter annihilation in the Sun and Earth respec- 
tively. In the case of the Earth, more than 12 events are 
needed for a DM mass between 70 GeV and 4 TeV. In 
the case of the Sun the number of events needed drops 
linearly as a function of mass starting from 300 events at 
70 GeV down to 70 events at 300 GeV. Beyond 300 GeV 
up to 4 TeV, the number of events needed remains fixed 
at 70. This provides a hint that one could detect the 
event rates caused by sneutrinos despite some differences 
between the sneutrino and neutralino dark matter spec- 
tra. These differences are due to unequal numbers and 
weighting of neutrino production channels, but the some- 
what harder spectrum of the sneutrino model will make 
IceCube slightly more sensitive to the model. Hence, we 
can expect that it might be possible to detect muon neu- 
trinos produced by sneutrino annihilation for sneutrino 
masses around 300 GeV for the Sun, cf. Fig. |4j Note that 
a large range of masses would be accessible with only an 
order of magnitude improvement in sensitivity. 

In summary, if the dark matter mass is determined 
from measurements at the LHC, then we can read the 
maximum number expected for the Sun muon rate in the 
B — L model from Fig. [If . Thus, for a known sneutrino 



Since we have used the upper bound on the sneutrino-proton 
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mass, observation of a muon signal exceeding the num- 
ber given in Fig. [4] will rule out the B — L model. The 
largest number of muon events from the Sun in the entire 
depicted mass range is 58 km~ 2 yr" 1 (36 km -2 yr _1 ) for 
case 1 (case 2). Therefore detection of a muon signal 
larger than this will rule out the B — L model regardless 
of the sneutrino mass. 

In the case of the Earth, as mentioned in the previous 
subsection, there is no prospect for a potential detection 
at IceCube for the standard halo model. The number of 
muon events is 6 orders of magnitude below the minimum 
measurable Earth rate of 12 km -2 yr" 1 in this case. 



V. DARK MATTER DISC IN THE MILKY WAY 

In our analysis, so far we have assumed a Gaussian like 
velocity distribution for dark matter particles with a typ- 
ical value for the three dimensional velocity dispersion of 
a v = 270 km sec -1 and |us un | = 220 km sec -1 for the ve- 
locity of the solar system with respect to the halo. How- 
ever, there are recent speculations about the existence of 
a dark matter thick disc in the Milky Way in addition to 
the baryonic one, see e. g. 0, Hf|. This dark matter disc 
is caused by the accretion of Milky Way satellite galax- 
ies and their corresponding baryonic and dark matter. 
As dynamical friction causes the satellite galaxies to ac- 
crete onto the disc, tidal forces disrupt the satellites [25| . 
Galaxy formation simulations find the density of the dark 
matter disc Pdark to be in the range ss 0.25—1.5 times the 
local halo dark matter density phaio [25[ .Possible ranges 
for the solar system velocity and velocity dispersion of 
the dark matter disc are|ug un | w — 150 km sec -1 and 
a v « 87 — 156 km sec -1 . 

Fig- HI shows the Earth muon rate when we scan about 
the relevant parameter space for the allowed values of 
I^Sunl and u v in case 1 and case 2. We used the fixed 
ratio Pdark/phaio = 1- Case 2 has a sufficient total event 
rate (> 12 km -2 yr _1 ) for nearly the whole allowed pa- 
rameter space. The constraint of the parameter space is 
more pronounced for case 1. The allowed combinations 
are roughly given by a triangle with maximal values of 
I^Sunl = 47km/s and a v — lOOkm/s. The differences in 
the allowed parameter space for the two cases reflects the 
fact that the Earth is not in equilibrium yet. Thus the 
muon neutrino signal and the corresponding muon flux 
still depends on the annihilation cross section, which is 
three orders of magnitude larger for case 2. However, we 
see that in both cases the Earth rates have increased to 
detectable rates, several orders of magnitude higher than 
the standard halo model that has higher |us«n| and a v , 
used in the previous section 11 . 
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FIG. 5: Total Earth-annihilation muon event rates inside the 
detector per kilometer squared per year for a 300 GeV (in 
case 1) and 1000 GeV (in case 2) sneutrino. 



A change in the velocities and dispersions also modifies 
the corresponding total Sun event rates. This is compar- 
atively modest for neutralino dark matter where it is at 
most one order of magnitude, see [26j]. Fig. [6] shows a 
band of allowed total Sun muon rates for the sneutrino 
dark matter. These rates are given again as a function 
of the sneutrino mass and under the requirement that we 



scattering cross section in our calculations, the number of muon 

events cannot be larger than that given in Fig. [4] 

The usage of a free space Gaussian velocity distribution means 



that our calculated event rates are an upper bound. There are 
many proposed parameterizations for the dark matter velocity 
distribution, and a Gaussian distribution belongs to the scenarios 
with the highest resulting event rates, see [2r3| . 
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have a measurable Earth rate of at least 12 events km - 
yr . Any variation of the event numbers for a fixed 
mass arises as a result of the use of velocities |i>s U n| and 
dispersions a v within the required parameter ranges of 
Fig. [5(a)] and [5(b)] A comparison between Fig. Q] and [5] 
shows an increase of ~ 30% in the Sun muon rates for 
the sneutrino dark matter. 




2000 
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FIG. 6: Total Sun-annihilation muon rates inside the detector 
for the sneutrino dark matter with modified velocity distribu- 
tions that yield Earth-annihilation rates of at least 12 events 
per year per km 2 . Upper (lower) curve shows the case 1 
(case 2). The dotted lines denote the mass range where one 
cannot explain the PAMELA data using case 2 anymore. 

The band of total muon rates for case 1 is notice- 



ably thinner than for case 2. It is seen from Figs. 5(a 



and 5(b) that case 1 has a smaller allowed parameter 



range with more than 12 events km -2 yr -1 . Thus the 
corresponding ratio between the minimal and maximal 
value within the allowed range is much smaller than that 
in case 2, and the possible change in the total Sun rates 
is comparatively small. Even for case 2 the differences 
between the highest and lowest rates for a fixed mass is 
about 40% or less. 

To summarize, a modified velocity distribution can 
substantially enhance the Earth muon rate for the sneu- 
trino dark matter beyond the detection threshold of 12 
km~ 2 yr -1 . It also raises the maximum Sun muon rate 
to 78 events km -2 yr -1 (48 km -2 yr -1 ) in case 1 (case 
2). Observation of the Sun muon rates larger than these 
will rule out the B — L model regardless of the sneutrino 
mass or Earth rates. 



VI. COMPARISON WITH MSUGRA 

Minimal supergravity (mSUGRA) is a constrained ver- 
sion of the MSSM that depends only on four parameters 
and one sign. These are mo (the universal soft breaking 
mass at the grand unification scale), mi/2 (the univer- 
sal gaugino soft breaking mass at the grand unification 
scale), A (the universal trilinear soft breaking mass at 
the grand unification scale), tan (3 (the ratio of MSSM 
Higgs VEVs at the electroweak scale) and the sign of /x 



(the MSSM Higgs mixing parameter). The mSUGRA 
dark matter candidate is the lightest neutralino. 

The parameter space of the mSUGRA model has 
three distinct regions allowed by the dark matter con- 
straints [291 ] : (i) the co-annihilation region where both mo 
and m.1/2 can be small, (ii) the hyperbolic branch/focus 
point region where the dark matter has a large Higgsino 
component and mo is very large but m 1 /2 is small, and 
(iii) the funnel region where both mo and m-1/2 are large 
and the dark matter annihilation occurs through heavy 
Higgs bosons in the s-channel. We note that a bulk re- 
gion (where none of the above properties hold) is now 
almost ruled out due to other experimental constraints. 
Among these three regions, the neutralino has a large 
capture rate in the hyperbolic branch/focus point re- 
gion due to a large Higgsino component that results in 
a large spin-dependent scattering cross section via Z ex- 
change. In this section we compare mSUGRA hyperbolic 
branch/focus point scenarios with the B — L model. 
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FIG. 7: Total Sun-annihilation muon rates inside the de- 
tector for mSUGRA hyperbolic branch/focus point sce- 
narios as a function of the neutralino mass. The results 
are for one year of detection with IceCube. 



Fig. shows the total Sun muon rate as a func- 
tion of the neutralino mass for mSUGRA hyperbolic 
branch/focus points. A comparison with Fig. [4] shows 
that these scenarios always have a higher total muon rate 
in the plotted mass range than the B — L model. The 
hyperbolic branch/focus point models yield larger muon 
rates by between more than one order of magnitude and a 
factor of 1.5 for dark matter masses in the 100 — 800 GeV 
range. Even for masses up to 400 GeV the hyperbolic 
branch/focus point scenarios provide rates higher than 
100 events km -2 yr -1 . These higher rates are explained 
by the bigger spin-dependent scattering cross sections, 
which are a few orders of magnitude larger than the up- 
per bound on the spin-independent cross section for the 
sneutrino dark matter. The spin-dependent scattering 
cross section for the B—L model is zero because U{\)b-l 
is a vectorial symmetry. Since the Sun mainly consists 
of hydrogen, the spin-dependent piece contributes domi- 
nantly for the mSUGRA case. 

However, it is interesting that despite having a much 
smaller scattering cross section, the B—L model can yield 
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muon rates that are roughly comparable to the mSUGRA 
scenarios. Sneutrino annihilation dominantly produces 
leptons, i.e., RH neutrinos in case 1 and taus in case 
2, which subsequently decay to LH neutrinos 100%. On 
the other hand, neutralino annihilation in the hyperbolic 
branch/focus point scenarios dominantly produces quark 
final states that have a small branching ratio for decay 
to neutrinos. 

Furthermore, despite lower event rates, sneutrino dark 
matter still produces a distinctive linear spectrum in the 
muon flux. As illustrated in subsection llVBl this feature 
is caused by the delta function in energy for the neutrino 
spectrum and can be used to distinguish between the B — 
L model and the hyperbolic branch/focus point scenarios 
as long as energy binning of the differential muon rate 
with respect to the energy is precise enough at IceCube. 
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FIG. 8: Total muon rates detected inside the Earth for 
mSUGRA Focus point scenarios as a function of the neu- 
tralino mass in the case of neutralino DM capture and 
annihilation in the Sun. Rates for a range of velocities 
and dispersions for which the corresponding Earth rates 
are at least 12 events per year per km 2 are shown in the 
shaded region. The results are for one year of detection 
with IceCube. 

Fig. [5] shows the counterpart of Fig. [5] for mSUGRA 
hyperbolic branch/focus point scenarios. The range of 
velocities and dispersions for which the corresponding 
Earth rates are at least 12 events km~ 2 yr _1 yields a 
band for the total Sun muon rates. We see that the 
range between the highest and lowest rates for a fixed 
mass does not exceed a factor of two even for masses 
below 200 GeV. 

A scan about the whole parameter space of the modi- 
fied velocity distribution yields a maximum of 13 events 
km~ yr _1 from the Sun for a 1000 GeV neutralino in 
the hyperbolic branch/focus point sneutrino. The B — L 
model with sneutrino masses of 1000 GeV and 1500 GeV 
gives rise to maximum values of 18 and 6 (25 and 14) 
events km -2 yr -1 for case 1 (case 2). In contrast, for a 
dark matter mass of 300 GeV, the maximum events km -2 
yr^ 1 are 158 (hyperbolic branch/focus point), 79 (case 
1) and 48 (case 2). Thus the hyperbolic branch/focus 
point rates are larger than the B — L rates for low masses, 
but both are in the detectable range at IceCube. At high 



masses it becomes more difficult to distinguish between 
the hyperbolic branch/focus point and the B — L models 
using maximal Sun rates; we would have to depend in- 
stead on the spectral features mentioned in Section llVBI 

In the stau co-annihilation and Higgs resonance re- 
gions the lightest neutralino has a high gaugino fraction 
and therefore a much smaller spin-dependent cross sec- 
tion that leads to much lower event rates than the B — L 
model. For example, even if we assume a modified ve- 
locity distribution without any minimal Earth event rate 
condition the maximum total Sun rate is less than 1 event 
km~ 2 yr^ 1 for a 300 GeV neutralino (compared with the 
maximum Sun rate of 158 events km -2 yr" 1 for a hyper- 
bolic branch/focus point scenario with the same mass). 
This is far below any detection threshold. 

It is also important to note that the hyperbolic 
branch/focus point in the mSUGRA model is incompati- 
ble with the g — 2 data, where there exists a 3<r deviation 
from the SM value if the e + e~ data is used to calculate 
the leading order hadronic contribution [3(| ■ In the con- 
text of the B — L model, case 2, which can address the 
PAMELA puzzle, also becomes incompatible with g — 2 
data, however the generic B — L model, i.e. case 1, is 
still compatible. 



VII. CONCLUSION 

We have considered prospects of indirect detection of 
the RH sneutrino dark matter in a U{1)b-l extension 
of the MSSM at the IceCube neutrino telescope. The 
sneutrinos captured in the Sun and Earth dominantly 
annihilate through S'-wave processes at the present time. 
In a generic situation (called case 1) the sneutrinos an- 
nihilate to RH neutrinos (annihilation cross section of 
3 x 10~ 26 cm 3 /sec) that quickly decay to a LH neutrino 
and the MSSM Higgs. If one seeks an explanation for 
the recently observed positron excess from the PAMELA 
data (called case 2), the sneutrinos with a mass > 1 
TeV dominantly annihilate to the lightest Higgs in the 
B — L sector (with an enhanced annihilation cross sec- 
tion of 3 x 10 -23 cm 3 /sec) that rapidly decay to fermion- 
antifermion pairs (74% taus, 16% bottom quarks, and 
10% RH neutrinos). LH neutrinos are produced mainly 
from the three-body decay of taus. The muon neutrinos 
from sneutrino annihilation are converted to muons via 
charged current interactions at IceCube. 

In both of the cases, sneutrino capture and annihi- 
lation inside the Sun reaches equilibrium. Consequently, 
the flux of neutrinos from the Sun is governed by the cross 
section for sneutrino-proton elastic scattering, which has 
an upper bound of 8 x 10~ 9 pb from the LEP and Teva- 
tron limits on the Z 1 mass (due to the vectorial nature of 
the B — L symmetry, there is no spin-dependent piece). 
In Fig. [4] we have shown the number of Sun muon events 
at IceCube as a function of the sneutrino mass for case 
1 and case 2 (using the upper bound on the sneutrino- 
proton scattering cross section). In both cases, the num- 
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ber of events are potentially detectable by IceCube due to 
a harder neutrino spectrum. Thus once the dark matter 
mass is found from measurements at the LHC, observa- 
tion of muon events larger than that given in Fig. [4] will 
rule out the B — L model. 

For the standard halo model the capture and anni- 
hilation of sneutrinos inside the Earth does not reach 
equilibrium for either case 1 or case 2, resulting in an 
event rate that is too small to be detected at IceCube. 
However, modified velocity distributions within the range 
allowed by recent simulations of the galaxy can lead to 
a substantially larger rate that exceeds the IceCube de- 
tection threshold of 12 km -2 yr _1 for events from anni- 
hilation in the Earth. Nevertheless, the Sun-annihilation 
muon rate can at most increase by 30% for a modified ve- 
locity distribution, as shown in Fig. [5] This implies that 
observation of a muon event rate larger than roughly 100 
km~ 2 yr _1 from the Sun will all but rule out the B — L 
model regardless of the dark matter mass. 

We compared predictions of the sneutrino dark mat- 
ter in the B — L model with that of the neutralino 
dark matter in the mSUGRA model. Only hyperbolic 
branch/focus point scenarios in mSUGRA, which have a 
Higgsino type dark matter candidate and thus large spin- 
dependent contributions to the neutralino-proton elastic 
scattering cross section, give rise to Sun muon event rates 
that can be detected at IceCube. Even though scatter- 
ing cross sections can be two to three orders of magnitude 
larger than the B — L case, the muon rates do not scale 
directly with the cross section. This is because sneutri- 
nos mainly annihilate into lepton final states (by virtue of 
the B — L symmetry) that decay to neutrinos with 100% 
efficiency, while neutralino annihilation dominantly pro- 
duces quark final states that have a small branching ratio 
for decay to neutrinos. Moreover, the linear dependence 
of the muon spectrum on the energy in the case of the 
sneutrino dark matter (particularly case 1) , a common 
feature for neutrinos produced from the two-body decays, 
can be used to distinguish between the B — L model and 
the hyperbolic branch/focus point scenarios. This will 
be feasible by a sufficiently precise energy binning of the 
differential muon rate at IceCube. 



APPENDIX A: MASS DEPENDENCE OF T A 

We analyze in detail here the contribution from Eq. ([3]) 
to the mass behavior of Fig. [4j Any mass dependence 
from dNl jdE v is ultimately washed out of the muon sig- 
nal by the linear dependence of (?cc,nc on the neutrino 
energy, which dominates at low energy. The distance D 
between the detector and the source and the branching 
fraction B~ into the final state / are independent of m^. 
Thus, the annihilation rate at high energies is governed 
by the dependence of r a on the kinematic suppression 
factor. Eq. ((3|) shows that this annihilation rate is pro- 
portional to the capture rate C. Ref. [I?) provides a 



parameterization of C as a function of the energy: 

C( TO w) « ^ F i( m N) S ( m N/ m N l ) lJ i( m N) > ( A1 ) 

i 

where the sum runs over all species i of nuclei in the 
Sun or Earth, Fi are the corresponding form factors, S 
is the kinematic suppression factor for capture of a sneu- 
trino and the <7j are the individual scalar cross sections 
for scattering from nucleus i. The effect of the Fi depen- 
dence on mass is negligible because most of these form 
factors vary little from unity. Furthermore, <7j is not de- 
pendent on the sneutrino mass in the B — L model since 
we have chosen a constant sneutrino-proton scattering 
cross section of 8 x 10 -9 pb (the upper bound implied by 
the Z 1 mass limits). Thus the overall shape of the curves 
in Fig. 2] can be understood by looking at S(m^). 
S can be parameterized by 



S(x) 



A(x) 



1.5 



i + A(xy- 5 



A(x) = 1.5 



{x- If 



3/2 



<V e , r >' 



(A2) 
(A3) 



where v — 270 km sec -1 is the velocity dispersion of the 
dark matter particles and < t> csc > is the escape velocity 
of 1156 km sec -1 and 13.2 km sec -1 for the Sun and 
Earth respectively. S(x) is bounded between zero and 
one. Moreover, it scales like 1.5(< v csc > 2 /v 2 )/x for 
x — > oo, and it peaks at one for x = 1. Therefore, the 
exact location of the peak for each scattering element 
i is determined by its corresponding nucleus mass m^r. 
As mentioned in subsection llV C[ S scales approximately 
like for large masses. The meaning of large in this 

context depends on the value of the (< v csc > /v) 2 ratio 
in comparison to x. For example, a value of with 
mjvf/mjv. > (< v csc> /v) 2 is considered large. 



APPENDIX B: CASCADE SIGNAL 

Fig. [HI plots the total energy spectrum from all cas- 
cades, both hadronic and electromagnetic (excluding the 
electromagnetic muon signal 12 ), per kilometer squared 
of detector per year for case 1 and case 2. The gen- 
eral downward trend of the plot occurs because the 
hadronic signal dominates as it is produced by both 
charged current and neutral current interactions, while 
the upward-trending electromagnetic signal only receives 
contributions from the charged current interactions and 
excludes the muon signal altogether. The cross sections 
for hadronic processes decrease as the transferred energy 



12 The electromagnetic cascade from a muon signal is excluded from 
the graph since it is accompanied by a more discernable muon 
track, the subject of the body of this paper. 
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FIG. 9: Total electromagnetic and hadronic cascades inside 
the detector volume from sneutrino annihilation in the Sun. 



to the nucleus goes up, hence creating the decreasing 
trend (high energy in the hadronic cascades corresponds 
to low energies in the electromagnetic cascades). We note 
that in the lower panel of the figure (case 2) the cascade 
signal is depleted at high energies. This is because the 
produced neutrinos have higher energies (as a result of 
the higher sneutrino mass in this case), and therefore ab- 
sorption and scattering effects inside the Sun are more 
important. This explains why the signal in case 2 is 
more steeply curved than the case 1 signal. 

It is important to remember that it is not clear at this 
time whether IceCube will be able to distinguish between 
electromagnetic and hadronic cascades. As a result, while 
a single charged current interaction will result in both 
a hadronic and electromagnetic cascade, these may be 
recorded as a single event with the total energy of the 
incoming neutrino. Meanwhile, the hadronic cascade of 
neutral current events would be recorded correctly as a 
single event with only part of the energy of the incoming 
neutrino. While we have assumed in the above that indi- 
vidual cascade signals are separable this may not reflect 
experimental reality. 
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